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A reduction procedure to obtain ground states of spin glasses on sparse graphs is developed 
and tested on the hierarchical lattice associated with the Migdal-Kadanoff approximation for low- 
dimensional lattices. While more generally applicable, these rules here lead to a complete reduction 
of the lattice. The stiffness exponent governing the scaling of the defect energy AE with system 
size L, a(AE) ~ L v , is obtained as yz = 0.25546(3) by reducing the equivalent of lattices up to 
L = 2 100 in d = 3, and as y 4 = 0.76382(4) for up to L = 2 35 in d = 4. The reduction rules allow 
the exact determination of the ground state energy, entropy, and also provide an approximation to 
the overlap distribution. With these methods, some well-know and some new features of diluted 
hierarchical lattices are calculated. 
PACS number(s): 05.50.+q, 75.10.Nr, 02.60. Pn. 



I. INTRODUCTION 

We propose a set of reduction rules applicable to spin 
glasses at T — on any sparse graph and arbitrary bond 
distribution. These reductions strip graphs of all vari- 
ables that are connected to at most two neighbors while 
accounting exactly for the ground state energy, entropy, 
and approximately for the overlap distribution of the sys- 
tem. In this paper, we introduce these reduction rules, 
and test them on the hierarchical lattice (see Fig. ^ ob- 
tained from the Migdal-Kadanoff bond-moving scheme 
0. The recursive nature of the hierarchical lattice per- 
mits us to quickly reduce the equivalent of 10 9 graphs 



with lengths corresponding to L 



>ioo 



d = 3 and 



L = 2 in d = 4 dimensions, limited only by accumu- 
lating rounding errors. Throughout this paper we focus 
exclusively on d = 3 and 4, and use only a discrete bond 
distribution, 



P(J)=pS(J 2 -l) + (l-p)5(J). 



(1) 



But our procedure is equally applicable to any continuous 
bond distribution, such as a Gaussian. 

The reduction produces high-accuracy results for the 
scaling of the defect energy width <j(AE) with L. We 
also study the diluted lattice, where p shall refer to the 
fraction of occupied bonds, and show that this scaling 
emerges for all p > p*. The critical point p* — 0.31032 2] 
is particular to the bond distribution in Eq. Q and lo- 
cated just above the percolation point p c = 0.281837 of 
the lattice. Since p c is a purely topological property of 
the lattice itself, having p > p c is merely a necessary con- 
dition for long-range correlated behavior in a spin glass. 
Cooperative effects from the bond disorder suppress cor- 
relations even for p c < p < p* , as was already discussed in 
Ref. |2j. Correspondingly, the moment (\q\) of the over- 
lap distribution |3j becomes non-zero only at p*, while 
the ground state energy and entropy remain smooth for 
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FIG. 1: Recursive generation of the hierarchical lattice, pro- 
ceeding from left to right. The reduction rules proceed from 
right to left, replacing the sub-graph with a reduced bond. 



all < p < 1, even at the transition. Both, entropy and 
overlap, exhibit an extremum at some p > p* which can 
be explained in terms of the peculiar lattice hierarchy. 

This study suggests that the reduction method pro- 
vides a powerful means to explore the (potential) on- 
set of replica symmetry breaking (RSB) 3] on finite- 
dimensional spin glasses for arbitrary bond distributions 
at T = but variable p. The onset of glassy behavior (al- 
beit not RSB) in the Migdal-Kadanoff model at p*, close 
to and intrinsically linked with p c , may suggest such a 
close link between percolation and the possible onset of 
RSB in realistic lattice models, which we are currently 
exploring £|. In fact, we already found a transition for 
the stiffness exponent in real three- and four-dimensional 
Edwards- Anderson models, resembling the one observed 
here, at a p* > p c 

Considered in terms of the average connectivity a = 
2dp ~ 0(1) of spins, dilute lattice glasses as well as the 
hierarchical lattice at T = behave strikingly similar to 
many mean- field models, which are currently studied in 
the context of combinatorial optimization ja, Q ■ Those 
models also enter into a glassy state (exhibiting RSB) at 
some finite connectivity a*, above a finite-connectivity 
percolation transition a c < a*. But in mean-field mod- 
els the onset of RSB is often studied as a precursor to 
yet another critical point, the satisfiability transition 0, 
whereas geometrically defined models arc unsatisfiable 
at any finite connectivity. These features should be in- 
dependent of the bond distribution, assuming zero mean 
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and unit variance. 

While the hierarchical lattice is completely reducible, 
diluted realistic lattices will become irreducible soon 
above percolation, in which case our algorithm || is used 
to reduce any graph as much as possible first, followed 
by a complete exploration of the ground states of the far 
more compact remainder graph using the extremal opti- 
mization heuristic 10]. With this hybrid approach, we 
are typically able to compute ground state energies, en- 
tropies and overlaps for graphs 10-100 times larger than 
previously recorded, extending well beyond the percola- 
tion point. 

II. MIGDAL-KADANOFF HIERARCHICAL 
LATTICE 

To illustrate the reduction method, and to test our 
algorithm, we consider here the hierarchical lattice 
(see Fig. QJ, obtained in the Migdal-Kadanoff real- 
space renormalization scheme []j for low-dimensional 
spin glasses. These lattices have a simple recursive, yet 
geometric, structure and are well-studied 
[3,0]. Most importantly, they are completely reducible, 
and we can discuss the method independently of any sub- 
sequent optimization that may be required for more com- 
plicated models. The most interesting property of these 
lattices is the curious fact that the scaling of its defect en- 
ergy distribution behaves very similar to that measured 
for actual two- and three-dimensional lattices [3, la, Il5l ] . 

As described in Fig. starting from generation 7 = 
with a single link, at each subsequent generation 7 + 1, 
all links from 7 are replace with a new subgraph. The 
structure of the subgraph arises from the bond-moving 
scheme in d dimensions, and has 2 d = 8 links for d = 3 
here. Thus, a hierarchical lattice of generation 7 has 
li = (2 d y links, thus corresponding to a d-dimensional 
lattice of "length" L = 2 1 but nj = 2+2 d - 1 (L d -l)/{2 d - 
1) = 0(L d ) vertices. While the average connectivity is 
2li/nj ~ 4 — 2 2 ~ d , the two root-vertices from genera- 
tion 7 = themselves obtain in generation J » 1 a 
connectivity of ~ 2( d ~ 1 ) 7 , and ~ 2 dI ~ 1 vertices are only 
two-connected, i. e. 7 in 8 for d = 3. 

The diluted hierarchical lattice percolates when there 
is a path between the two root- vertices. This notion leads 
to a simple recursion relation for the percolation thresh- 
old by counting the weights of all diluted subgraphs from 
Fig. [I] that percolate: 

p I+1 =4pj-6pj + Ap 6 I -pl (2) 

which has a non-trivial stationary point at _p c = 
0.2818376366. It has been pointed out by Ref. that 
a spin glass with the discrete ±J-bond distribution in 
Eq. Q exhibits instead a critical transition between a 
paramagnetic and a spin glass phase at p* = 0.31032, 
which is closely related to the percolation transition. 
While below p c disconnected bonds (J = 0) clearly 
dominate and prevent long-range correlations, even for 



p c < p < p* such correlations remain suppressed due to 
the cooperative behavior of parallel bond structure per- 
vasive in the lattice that lead to many cancellations (see 
"double rule" in Sec. II I II) and additionally disconnects 
subgraphs at some higher level of the hierarchy. 

III. REDUCTION RULES FOR SPIN GLASSES 

In the Ising spin glass problem we assign to each ver- 
tex i of a graph or lattice, 1 < i < n, a spin variable 
xt e {— 1, +1}, while each link between two connected 
vertices i and j obtains a bond variable Jij that is drawn 
at random from some distribution P(J), which may be 
discrete or continuous. The problem consists of find- 
ing spin configurations {x\, . . . ,x n ) which minimize the 
Hamiltonian, or total energy, 

H{x\ . . ■ , 37 n ) — ^ 7i,j xi Xj (3) 
<ij> 

for a fixed (quenched) set of bonds Ji.j. Of course, for 
Ji.j > (only ferromagnetic bonds), the solution is sim- 
ply to make all connected spins aligned with each other. 
But if some or all bonds are negative, and spins are suffi- 
ciently connected, the problem can be frustrated 0] and 
solutions become nontrivial. 

In the following, we explain the rules by which to 
eliminate recursively all two-connected spins while ac- 
counting for the energy, entropy, and approximate over- 
lap of ground-state configurations. To this end we have 
to generalize bonds by adding internal degrees of free- 
dom, which evolve during the reduction process. Each 
bond between two spins i and j now consists of a tuple 

^id = \ ' m i,j ' m i ,j ' S i ,j ' S i,j) ' W 

where Jj j is the weight of the bond, m i j is the entropy 
accumulated by the bond, and s i • counts the number of 
previously reduced spins whose state is completely deter- 
mined ("entrained") by the bond for XiXj = ±1, respec- 
tively. In addition, there is an energy offset -Eoffsct _ 
accounting for the energy difference between the original 
and the reduced graph. In general, determining entropy 
and overlap from the m- and s-values can be quite com- 
plicated for the interesting case of an irreducible remain- 
der graph 0. For the completely reducible hierarchical 
lattice, it is a lot simpler and we will describe it below. 
Initially, for the unreduced graph, all m- and s-values 
and -Eoffsct are zero, and the bonds are drawn from one 
of the usual distributions, for example, Eq. l|T]l. 

These rules are elaborate but easily enumerated, inde- 
pendent of the previous or future structure of the graph 
or the spin configuration. A more complete set of rules 
including zero-, one-, and three-connected vertices re- 
quired for arbitrary graphs will be presented elsewhere 
0, There, each spin i obtains one additional num- 
ber Si (aside from its state Xi) that counts the previously 
eliminated spins whose state is entrained to i. 
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To demonstrate the use of these extra degrees of free- 
dom, consider, for instance, the reduction of a one- 
connected spin j: It and its bond J\ j to its sole neigh- 
bor i disappears from the graph, and the connectivity 
of i is decremented. The new graph will have itself a 
ground state energy offset by \Jij\ relative to the unre- 
duced graph, i. e. -Eoffset <— -^offset - \Ji,j\- For J i:j = 
(a possibility we eventually have to allow for, even if 
our unreduced bond distribution did not), all information 
contained in Jij is discarded and s, doesn't change, since 
Xj is free to take on any value ±1 in the ground state; 
only the entropy of the old graph is offset by In 2 with 
respect to the new graph. In turn, for Jjj > ( Jij < 0), 
j is entrained to the state of i in any ground state con- 
figuration, and Si is incremented. If the bond Ji t j had 
previously inherited nonzero entropies and entrain- 
ments •, there would be an offset in entropy by mfj 
(m~ ) between the reduced and the original graph, and Si 



would be further increased by s 



+ 



m 



+ 



) and 



1,3 \ *,3'i 1,3 V %3 > 

s 7j ( s tj) are discarded because the ground state must 
have XiXj = +1 (xtXj = — 1). 

With the two-point rule we use exclusively in this pa- 
per, a spin variable Xi and its two bonds J\ and Ji are 
replaced by a new bond J . It is easy to show that 



J 



(5) 



which assigns to Jj^k the bond with the smaller modulus 
and the sign of J\ J 2 . The energy offset E Q g set , account- 
ing for the energy of the original graph, will have to be 
lowered by 



\{\Jx + M + \Ji-M 



(0) 



At this level, our reduction rules correspond to the fa- 
miliar traces (at T = 0) used many times before to study 
the hierarchical lattice (see, e. g., Refs. 0,0). 

Clearly, for | Ji| ^ | J 2 |, the reduced spin i is entrained 
to the neighbor with the bond of larger modulus, while 
the weaker bond may be violated. The entropy and en- 
trainment variables that J inherits from the two reduced 
bonds arc controlled by the sign of the stronger one. We 
can summarize these rules by defining J7> and J < as the 
bond with the stronger and weaker (in modulus) weight, 
respectively. Then, we obtain for the new bond: 



\J>\ > \J< 



4- siqn(J > ) . ±sinn(J>) 



sign(J > ) ±sign(J > ) 



The case \J\ \ = J 2 | leads to further distinctions: 



(7) 



sign(J!) sign(.h) 



Jl = ±J 2 + 




^sign(Ji) sign{Jo.) 



, Ji = ± J 2 ^ 



(8) 



Finally, we need a rule to combine double bonds, which 
will be created eventually. When two bonds, J\ and J%, 
connect the same two spins, they are simply replaced 
with a new bond 



J = J1+J2. 



(9) 



It is important to note that the two-point and double 
rules only effect the newly created bond and the energy 
offset, but not the two neighboring spins (or the entropy 
offset). Unlike for the one-point rule, we do not need to 
consider an internal entrainment number for spins here. 



IV. ALGORITHM AND IMPLEMENTATION 

With the reduction rules from Sec. IIHI it is now 

straightforward to evaluate hierarchical lattices of any 
size. In effect, these rules allow us to reverse the generat- 
ing mechanism depicted in Fig. ^ Recursively, from one 
generation to the next, we reduce all two-connected spins 
in the middle of the subgraph first, then use the double- 
bond rule sequentially, until a single, reduced bond is 
left. Once an entire graph has been reduced into such a 
single bond, the weight J of that bond is in fact half the 
defect energy, AE = 2 J, of that graph, i. e. the energy 
difference between having both root spins aligned and 
anti-aligned. The entropy of that graph is simply given 
by m sl9n ( J \ or by ln[exp(m + ) + exp(m - )] for J = 0. 
The energy of the graph can be obtained from the run- 
ning energy-offset -Eoffset , accounting for the reduction of 
bonds along the way, and J. For our purposes, it is more 
useful to extract the (non-negative) "cost" of a graph 
given by the absolute weight of all violated bonds for 
any bond distribution. In case of the discrete ± J dis- 
tribution in Eq. Q), this sum reduces to a count of all 
violate bonds (all of absolute weight 1) (l+E g set — \J\)/2, 
where (I) = pL d is the average number of non-vanishing 
bonds in a graph. As an approximate measure of the 
overlap, we merely store the largest entrainment number 
Smax observed at any point during the reduction process, 
which is a measure of the largest correlated cluster of 
spins within the lattice. The size of this cluster will ei- 
ther dominate the overlap (with smaller clusters adding 
some negligible fluctuations) in the spin glass state or be 
itself a vanishing fluctuation in the paramagnetic state. 

Although these measurements are not difficult, the re- 
duction of ever larger graphs soon becomes untenable, 
since the cost of reducing true graphs would grow by a 
factor of 2 d for each generation. Instead, we exploit the 
fact that each subgraph within the hierarchy is indepen- 
dent of any other parallel subgraph. As above, we reduce 
all graphs down to a single bond. To assemble a graph 
of generation 7 + 1 we only need a large enough pool 
of independent graphs of generation 7. In fact, many 
nearly-independent graphs of generation 7+1 can be as- 
sembled by repeatedly drawing at random on 2 d graphs 
of generation 7 from that pool. 
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We proceed as follows: Assume that we already possess 
a pool of Aj reduced graphs at every generation I (up 
to a maximal generation / max ), each graph represented 
by its bond J , E g sct , and s max . Each cycle of the al- 
gorithm generates an elementary subgraph (generation 
1=1) with bonds J drawn from the bond distribution 
in Eq. Q and empty entropy, entrainment, E ff sct , and 
Smax variables. The elementary subgraph is reduced to 
obtain a bond J of the first generation, replacing the 
oldest entry in pool A\. After every fcth new reduced 
bond entering into pool Ax—i, we randomly assemble 2 d 
such bonds into a new subgraph to create a new reduced 
bond to be entered into pool Ai, as long as I < / max - 

If we choose Aj = k = 2 d and assemble bonds se- 
quentially instead of at random, this algorithm produces 
exact hierarchical lattices, where each subgraph is certain 
to be independent. But to obtain just one graph of, say, 
generation / = 10, we have to assemble 2 30 elementary 
subgraphs for d — 3. Instead, for all p not too close to p*, 
we simply choose large pools, Aj = 2048, and as small 
as k = 1. That is, each time we assemble an elementary 
graphs at I = 0, we enter its value to pool A\, assem- 
ble a subgraph out of A\, add the reduced bond to pool 
A2, etc, all the way up to / max in each cycle. Hence, 
our implementation becomes independent of L, and for 
instance at p = 1 we have assembled graphs up to gen- 
eration I m ax = 100 in d = 3 and J max = 35 in d = 4, 
corresponding to L = 2 100 or 2 3 5, respectively. / max is 
only limited by apparent limits on floating point accu- 
racy due to repeated cancellation between almost equal 
bonds. Note that each pool gets completely refreshed 
after Ai creations of new graphs from the pool just be- 
low it, while there are about 0(2 dAl ) different subgraphs 
that could be created out of each pool at any one time. 

Every time a graph gets added to a pool, its proper- 
ties are evaluated from its J, E s set , and s ma x- While 
ground state energy, entropy, and (approximate) overlap 
converge quickly with J, and throughout are quite in- 
dependent of the choice of k and Aj 1 the measurement 
of the defect energy becomes quite sensitive to the algo- 
rithmic parameters near p* . In those cases, we increase 
k (and lower I max correspondingly), until stable results 
were obtained. Only at p w p* , the exact algorithm was 
necessary to observe the true defect energy distribution. 



V. NUMERICAL RESULTS 

Using the algorithm from Sec. II VI we have studied vari- 
ous properties of a Ising spin glass with bond distribution 
as given in Eq. JQ| on the diluted and undiluted hierar- 
chical lattice. In fact, using the easiest implementation 
with k = 1, pool sizes of Aj = 2048, and / max = 20, we 
first scanned lattices for < p < 1 in steps of Ap = 0.1. 
At each p and I, we averaged over 10 7 values. We found 
that for each p, the ground state energies per spin, en- 
tropies per spin, and the moment of the overlap, (\q\) 
were well-converged at J max = 20, and insensitive to the 
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FIG. 2: Plot of the width a(AE) of the defect energy distri- 
bution as a function of systems size L in d — 3 for various 
bond fractions p. Dashed lines from bottom to top refer to 
p = 0.1, 0.2, . . . , 1, while data with diamonds and crosses re- 
fer to p = 0.28, 0.29, . . . , 0.34. The data points with circles 
were obtained using the k = 1 implementation from Sec. II V I 
those with diamonds using k = 2, and those with crosses at 
p = 0.31 ~ p* required the exact algorithm. Note that for 
p < p* the width evolves toward the p = fix-point with 
<j(AE) = 0, while for p > p* it invariably evolves to the p = 1 
fix-point with scaling behavior cr(AE) ~ L v . While the exact 
asymptotic behavior emerges only for very large sizes L when 
p ~ p* , it is already obvious for small L whether a(AE) is 
going to rise or to fall. 



choice of k using trials with k = 2. 

We also found t he width of the d efect energy distri- 
bution, a(AE) = ^{AE 2 ) - (AE) 2 , to be stable except 
for p = 0.3. For large I = \og 2 (L) the distribution of 
defect energies either collapses into a ^-function at zero 
for small p, or eventually converges to the distribution of 
the undiluted lattice (p = 1). Near p w p* , the evolu- 
tion of that distribution becomes highly sensitive to small 
fluctuations, and the minute statistical dependencies in- 
herent in our algorithm for small k escalate. Thus, we 
ran another scan for 0.28 < p < 0.34 using k = 2 and 
10 4 graphs at I max = 20 (and hence gaining a factor k 
more graphs at each lower generation /). Now, only for 
p = 0.31 » p* did <t(AE) become unstable (by compar- 
ison with trial runs at k = 4). Finally, for p = 0.31, 
we used an exact algorithm and indeed found a(AE) 
to be horizontal for many generations until it eventually 
decays, indicating that p* > 0.31, consistent with the 
expected value of 0.31032 from Ref. 0- The results for 
a(AE) as a function of / = \og 2 (L) for all p are shown 
in Fig. 

In Figs. 13151 we plot the number of violated bonds 
(or cost) per spin, entropy per spin, and a characteris- 
tic moment of the (approximate) overlap distribution for 
L — > co and d = 3. The cost is already non-zero for small 
p, rising consistent with p 4 , proportional to the proba- 
bility of creating a small, frustrated loop that requires 
four bonds. Just beyond p* , the cost per spin rises more 
slowly, but without drastic changes for any p. Similarly, 
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FIG. 3: Log-log Plot of the cost per spin for the infinite system 
in d — 3 as a function of p. This cost, i. e. the fraction of 
violated bonds, only approaches zero for p — > (consistent 
with p 4 ). The data smoothly continues through both critical 
points at p c and p* , indicated by dashed vertical lines. For 
p > p* the cost appears to rise less rapidly than for p < p* . 



the entropy density in Fig. 0] remains unaffected by the 
transition at p* , as is expected Instead, it reaches a 
minimum at about p = 0.7 that is an artifact of the hier- 
archical lattice. The large majority of spins in the graph 
(a fraction ~ 7/8 for d = 3) are at most two-connected 
spins whose contribution dominates the entropy density. 
In a fully bonded graph (p = 1), about half of these spins 
are free, the other half is constraint, with an energy den- 
sity of ln(2)/2. For small dilutions p < 1, about every 
removed bond turns a two-connected spin into a one- 
connected spin, which is always constrained, and the en- 
tropy declines. Eventually, at large dilution (small p), an 
increasing number of one-connected spins become totally 
disconnected and free again, raising to entropy density 
to its ultimate value, ln(2), for the unconnected (p = 0) 
lattice. By this argument the entropy density should be 
given simply by 



S 
n 



ln(2), 



(10) 



which fits the data exceedingly well. 

In contrast with energy and entropy, the overlap distri- 
bution is clearly sensitive to the transition from paramag- 
netic to glassy behavior at p* , as shown in Fig.JSJ There, 
we plot only its moment, since the distribution itself 
quickly converges to a structureless (5-peak as expected 
for a model that is replica symmetric (RS) throughout. 
In fact, we find that (\q\) / •*/ (q 2 ) ~ 1 to within 1% for 
all p. (For one, a broad overlap distribution indicative of 
widely separated ground state configurations and RSB 3] 
would be antithetical to the reducibility of the hierar- 
chical lattice.) About at p*, the location of the peak 
becomes non-zero and splits into two symmetric peaks, 
leading to a non-zero moment (\q\). 

The moment experiences a surprising maximum very 
close to p = 1, varying sharply for small dilutions. 
Again, the maximum may be explained in terms of the 




FIG. 4: Plot of the entropy per spin for the infinite system as 
a function of p. The data smoothly continues through both 
critical points at p c and p* , indicated by dashed vertical lines, 
but exhibits a surprising minimum. The data is approximated 
well by Eq. 1101 (continuous line), as discussed in the text. 




FIG. 5: Plot of the moment (|g|) of the overlap distribution 
as a function of p. The moment remains zero for all p < p* , 
including at the percolation point p c . (Both, p c and p* , axe 
marked by dashed vertical lines.) Above the transition at p* 
to a spin glass state, (|q|) rises first linearly, but then turns 
over at a maximum at small dilutions just below p — 1, and 
drops sharply at p = 1. The inset shows that (|g|) seems to 
scale with (1 — p) 1//5 for p — *• 1; the dashed line is a linear fit. 



abundance of two-connected spins: Whatever the size of 
strongly correlated spin clusters may be in the fully con- 
nected lattice (p = 1), the overwhelming conversion of 
two-connected to one-connected spins at small dilution 
can only increase the size of those clusters at its edge, 
since one-connected spins are always entrained. But more 
importantly, the elimination of two-connected spins, say, 
within a subgraph as shown in Fig. ^ wm enhance the 
probability that the two root-spins of that subgraph be- 
come correlated, which would merge the clusters that 
those spins are part of. While the latter effect could 
suddenly increase the size of typical clusters by a factor, 
it is not clear whether the steep slope observed at p = 1 
is not an artifact of our approximation to the overlap (see 
Sec. IIV|I . Further numerical studies suggest that 



(M)p-<M)p=i~(i-p)« (p-»i), 



(ii) 
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as seen in the inset of Fig. [3] We have not found an 
explanation for the scaling in Eq. (|llf> in terms of the 
lattice yet, since it is also a property of the entrainment 
in the spin glass state (and hopefully, the overlap). 

As shown in Fig. |2 for all p > p* we obtain the ex- 
pected scaling of the width of the defect energy distribu- 
tion with system size L, 



a(AE) ~atf [1 + e(L)} 



(12) 



for large enough L, including a next-to-leading correc- 
tion e(L) <C 1. Most conveniently, it appears that our 
k = 1 algorithm (for large enough pools Aj) described 
in Sec. IIVI is least likely to create spurious dependen- 
cies for p = 1, where we would expect the quickest and 
most stable convergence toward the renormalized bond 
distribution. Thus, we disabled the sampling of energy- 
offsets, entropies, and entrainments, which already are 
well-converged at J max = 20, and only focused on the 
reduction of bonds. Setting 7 max = 100 (d — 3) or 
/ max = 40 (d = 4) , k = 1, and A/ = 1024 for each 
/, we measured more than 10 9 reduced bonds at each I 
to obtain the distributions of defect energies AE — 2 J 
in d = 3 and 4. (In fact, we repeated the same calcu- 
lation with a different random number generator, with- 
out markable difference in results. A previous run with 
I,„ax = 150 in d = 3 showed divergence due to limited 
floating point precision at 7 w 120; the same problem 
emerged already at I max = 40 in d — 4.) At that level, 
we obtained a ratio between its absolute mean and its 
width a(AE) of < 10 -4 in both dimensions, even for the 
largest /, which is a measure of the statistical error (the 
distribution is supposed to be symmetric). 

Surprisingly, although we extended the scaling regime 
in d = 3 by a power of 10 compared to previous calcula- 
tions (L = 2 100 compared to L = 2 10 [l3|), the extrapo- 
lation for the scaling exponent y in Eq. (|12|) only yields 
5 digits accuracy, see top of Fig. H3 First, we approxi- 
mated y « 0.2555 from the intercept of log 2 (er)/ log 2 (L) 
with respect to 1/ log 2 (L) at L — > oo, which displays large 
corrections to scaling. We are then able to refine our es- 
timate of y by converting Eq. (|12J) to ln(er) — y \n(L) ~ 
ln(a) + e(i), using ln(l + e) ~ e. Plotting that expres- 
sion again vs. l/log 2 (L) (note that log 2 (L) = /) sug- 
gest a straight-line asymptotic behavior that would im- 
ply e ~ bj log 2 (L) for L — > oo. The top of Fig.EJindicates 
that a value of y — 0.25546(3) in d — 3 leads to an appar- 
ently converging result with ln(a) « 3/4 at the intercept 
for L — > oo. We have used the same procedure to find 
y = 0.76382(4) in d = 4, see bottom of Fig.© Although 
the floating point accuracy here deteriorated already at 
L = 2 35 , the asymptotic scaling regime is reached for 
much small L than in d = 3, leading to comparable ac- 
curacy in the fit. 

Our value for the stiffness exponent y in d = 3 is 6% 
below y = [d-l+log 2 (l-2/7r)]/2 = 0.2698 obtained ana- 
lytically in Ref. |l2j | with a Gaussian approximation, and 
consistent with y « 0.25 obtained numcr ically in Ref. @ 
by evolving a "pool" of bonds through several genera- 
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FIG. 6: Extrapolation plot for the leading asymptotic be- 
havior of a(AE) as a function of system size L in d = 3 
(top) and d — 4 (bottom). The data for a(AE) is converted 
to ln(a) — yln(L) and plotted vs. the apparent scaling cor- 
rection, l/log 2 (I/), for different estimates of the exponent y. 
For d = 3 the graph at the center (circles) plots the data 
for y — 0.25546 and seems to converge linearly to a finite 
value of about 3/4 for L —* oo. In turn, using y = 0.25549 
(down triangle) or y — 0.25543 (up triangle) on the same 
data seems to lead tor diverging extrapolations, suggesting 
yz = 0.25546(3). The same procedure for d — 4 brackets the 
"best" choice of y — 0.76382 (circles) between two diverging 
choices of y — 0.76386 (down triangle) and y = 0.76378 (up 
triangle), suggesting 2/4 = 0.76382(4). 



tions, very similar to our algorithm in Sec. If VI Since 
both of these references use an initial Gaussian bond 
distribution, in contrast to our discrete distribution in 
Eq. Jp), the results suggest that the stiffness exponent is 
independent of the bond distribution. In Ref. [lj|i an ex- 
pansion in the limit of large dimension d — > 00 is used to 
demonstrate the independence of y of the initial bond dis- 
tribution. A recursion in the generation / is constructed 
which for / — > 00 leads invariably to the Gaussian de- 
rived in Ref. 01 ■ With finite-dimensional corrections, 
this expansion improves on the Gaussian result above to 
yield 2/3 ~ 0.255, in excellent agreement with our result 
here. In contrast, we expect that energy, entropy, over- 
lap, and even p* would depend on the distribution at 
least quantitatively. 
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VI. CONCLUSIONS 

We have presented a new method to reduce spin glasses 
on arbitrary graphical structures. The reductions are 
most appropriate for sparse graphs near their percola- 
tion threshold, which may be either completely (as in 
this paper) or substantially reduced, facilitating further 
optimization using standard methods [4J, |5j. In either 
case, the reductions allow to extract exact ground state 
energies and entropies, and even provide a good approx- 
imation to collective properties such as the overlap dis- 
tribution. 

With the reduction method we have reproduced a num- 
ber of ground-state properties of the Migdal-Kadanoff hi- 
erarchical lattice at any dilution, and discovered various 
new features, which can be mostly explained in terms of 
the lattice structure. The simplicity of this model further 
permitted us to determine the exponent y in d = 3 and 4 
for the growth of the defect energy with system size to 
unprecedented accuracy. We have recently applied the 



reduction method to diluted Edwards- Anderson models 
for d = 3 and 4, obtaining y 3 = 0.240(5]and y 4 = 0.60(1) 
for the respective stiffness exponents 5]. These results 
demonstrate that the Migdal-Kadanoff approximation in 
d = 3 is even closer to the Edwards- Anderson value than 
had been previously expected 0, 0] , although they are 
definitely distinct. The differing values in d = 4 reflect 
the fact that Migdal-Kadanoff is a low-dimensional ap- 
proximation. 

The results obtained here and in Ref. p| suggest that 
the reduction method may be a useful tool in the explo- 
ration of Ising spin glass systems on any graph or lattice 
near percolation that may reveal the onset of non-trivial 
configurations of ground states. 
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